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TRANSPORT  PHENOMENA  AND  INTERFACIAL  KINETICS 
IN  MULTIPHASE  COMBUSTION  SYSTEMS 


Principal  Investigator:  Prof.  Daniel  E.  Rosner 
1.  INTRODUCTION 

The  performance  of  ramjets  burning  slurry  fuels  (leading  to  condensed  oxide  aerosols 
and  liquid  film  deposits),  gas  turbine  engines  in  dusty  or  marine  atmospheres,  or  when  using 
fuels  from  non-traditional  sources,  depends  upon  the  formation  and  transport  of  small 
particles  across  non-isothermal  combustion  gas  boundary  layers  (BLs).  Even  airbreathing 
engines  burning  "clean"  hydrocarbon  fuels  can  experience  soot  formation/deposition 
problems  (e.g.,  combustor  liner  burnout,  accelerated  turbine  blade  erosion  and  "hot" 
corrosion).  Moreover,  particle  formation  and  transport  are  important  in  many  chemical 
reactors  used  to  synthesize  or  process  aerospace  materials  (turbine  blade  coatings,  optical 
waveguides,  ceramic  precursor  powders,  fibers  for  composites,...).  Accordingly,  our 
research  is  directed  toward  providing  chemical  propulsion  systems  engineers  and  materials- 
oriented  engineers  with  new  techniques  and  quantitative  information  on  important  particle- 
and  vapor-mass  transport  mechanisms  and  rates. 

The  purpose  of  this  report  is  to  summarize  our  research  methods  and 
accomplishments  under  AFOSR  Grant  91-0170  (Technical  Monitor:  J.M.  Tishkoff)  during 
the  3-year  period:  15  February  '91-14  February  '94.  Readers  interested  in  greater  detail 
than  contained  in  Section  2  are  advised  to  consult  the  published  papers  explicitly  cited  in 
Sections  2  and  5.  Copies  of  any  of  these  published  papers  (Section  5.2  and  Appendix)  or 
preprints  (Section  5.3)  can  be  obtained  by  writing  to  the  PI:  Prof.  Daniel  E.  Rosner,  at  the 
Department  of  Chemical  Engineering,  Yale  University,  New  Haven  ,CT  06520-8286  USA. 
Comments  on,  or  examples  of,  the  applications  of  our  research  (Section  3.4)  will  be 
especially  welcome. 

An  interactive  experimental/theoretical  approach  has  been  used  to  gain 
understanding  of  performance-limiting  chemical-,  and  mass/energy  transfer-phenomena  at 
or  near  interfaces.  This  included  the  development  and  exploitation  of  seeded  laboratory 
burners  (Section  2.1),  new  optical  diagnostic  techniques  (Section  2.2)  and  flow  reactors 
(Section  2.4).  Resulting  experimental  rate  data  ,  together  with  the  predictions  of  asymptotic 
theories  (Section  2),  were  used  as  the  basis  for  proposing  and  verifying  simple  viewpoints 
and  effective  engineering  correlations  wirth  a  rational  basis  for  future  design/optimization 
studies. 


2.  RESEARCH  ACCOMPLISHMENTS 

Most  of  the  results  we  have  obtained  under  Grant  AFOSR  91-0170  during  '91 -'94 
can  be  divided  into  the  subsections  below: 

2.1.  Transport  and  Stability  of  Aggregated  Particles:  Theory 

The  ability  to  reliably  predict  the  transport  properties  and  stability  of  aggregated 
flame-generated  particles  (carbonaceous  soot,  AI2O3,  Si02,...)  is  important  to  many 
technologies,  including  chemical  propulsion  and  refractory  materials  fabrication.  The 
existence  and  character  of  such  particles  is  also  known  to  influence  the  "signature"  of 
chemical  propulsion  devices. 


The  Brownian  diffusion-,  inertial-,  and  optical-properties  of  aggregated  particles,  as 
formed  in  sooting  diffusion  flames,  are  quite  sensitive  to  size  (e.g.  number  N  of  "primary” 
particles;  see  Fig.  1)  and  morphology  (geometrical  arrangement  of  the  primary  particles).  In 
this  program  we  developed  methods  to  anticipate  coagulation  and  deposition  rates  of 
suspended  populations  of  such  particlesjin  combustion  systems.  As  one  example,  we  have 
recently  developed  improved  and  efficient  methods  for  predicting  the  Stokes  drag  of  large 
’fractal'  aggregates  via  a  spatially  variable  porous  sphere  model  (Tandon  and  Rosner,1994; 
Figs.  1,2).  Using  the  Stokes-Einstein  equation,  the  results  of  Fig.  2  have  been  used  to 
predict  the  Brownian  diffusivity  of  such  aggregates  in  the  high  pressure  (near  continuum-) 
limit  (proportional  to  the  product  of  the  reciprocal  of  the  ordinate  of  Fig.  2  and  N-1/Df).  This 
approach  can  be  extended  to  predict  the  thermophoretic  diffusivity  of  such  aggregates,  an 
important  quantity  we  have  recently  found  to  be  much  less  sensitive  to  size  and  morphology 
than  the  translational  Brownian  diffusivity  (Rosner  et.al.  1992).  Indeed,  this  provides  the 
theoretical  basis  for  the  thermophoretic  sampling  technique  being  employed  in  our  current 
experimental  studies  (Section  2.2).  These  new  methods/results,  together  with  recent  results 
on  the  spread  of  aggregate  sizes  in  coagulating  populations,  can  be  used  to  predict  wall 
capture  rates  by  the  mechanisms  of  convective-diffusion,  turbulent  eddy-impaction,  and 
thermophoresis.  Also  developed  in  this  program  were  efficient  pseudo-continuum  methods 
to  predict  chemical  interactions  between  aggregates  and  their  surrounding  vapor 
environment— interactions  which  can  lead  to  primary  particle  growth,  or  burn-out.  In 
particular,  we  developed  new  and  efficient  methods  to  predict  the  "accessible  surface  area" 
of  aggregates  (expressed  as  a  fraction  ,  r\,  of  the  true  surface  area  in  Fig.  3),  including  its 
dependence  on  size  (N),  structure  (fractal  dimension,  Df),  probing  molecule  reaction 
probability  cl,  and  pressure  level  (via  Knudsen  number  based  on  primary  particle 
diameter)(Rosner  and  Tandon,  1994). 

Initiated  in  this  program  were  studies  of  the  restructuring  kinetics  of  aggregates  — 
te.  those  factors  which  determine  the  observed  size  of  the  apparent  "primary  particles" 
comprising  soot  particles,  and  the  "collapse"  of  surface  area  observed  in  some  high 
temperature  systems  (Cohen  and  Rosner,  1993).  Toward  this  end,  we  developed  new 
methods  to  characterize  the  morphology  of  multi-particle  aggregates  (Fig.  5) 
thermophoretically  extracted  from  laminar  CDFs  in  a  new  "slot"  type  burner  (Fig.  4).  One 
such  "fingerprint"  is  the  pdf  of  angles  formed  by  triplets  of  primary  particles  (Fig.  6  ). 

2.2.  Formation,  transport  and  stability  of  combustion-generated  particles- 

LAMINAR  COUNTERFLOW  DIFFUSION  FLAME  EXPERIMENTS 
(  We  have  inferred  the  thermophoretic  diffusivity  of  flame-generated  submicron 
"soot"  particles  using  two-phase  flame  structure  measurements  on  (TiCl4(g)-)seeded  low 
strain-rate  counterflow  laminar  diffusion  flames  (CDF-)  (Gomez  and  Rosner,  1993).  A 
knowledge  of  the  relative  positions  of  the  gas  and  particle  stagnation  planes  and  the 
associated  thermal  and  chemical  environments  can  be  used  to  control  the  composition  and 
morphology  of  flame-synthesized  particles.  These  factors  should  also  influence  particle 
production  and  radiation  from  turbulent  non-premixed  "sooting"  flames,  as  discussed  further 
in  Gomez  and  Rosner,  1993. 

To  obtain  fundamental  information  on  nucleation,  growth  and  aggregate  restructuring, 
we  developed  an  improved  "slot-type"  burner  (Fig.  4  )  and  introduced  instruments  to  carry 
out  in  situ  measurements  of  particle  Brownian  motion  ( via  "dynamic  light  scattering").  We 
also  developed  a  thermophoretic  sampler  to  extract  aggregates  from  various  positions  in  the 
seeded-CDF  for  morphological  analysis  using  transmission  electron  microscope  (TEM) 
images  (Fig.  5).  Aggregate  data  obtained  from  CH4  flames  seeded  with  titanium  tetra- 

isopropoxide  (TTIP-)  vapor  are  being  being  analyzed  using  new  theoretical  methods  briefly 
outlined  in  Sections  2.1,  2.3. 
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Fig.l  "Porous  sphere"  model  of  large  fractal  aggregate 
suspended  in  a  background  gas;  basis  for  the  calculation  of 
translational  an(j  rotational  Brownian  diffusion 
coefficients,  thermophoretic  diffusivity,  "stopping  time", 
accessible  area,  and  restructuring  kinetics  (after  Rosner 
and  Tandon,  1993,  Rosner,  Cohen  and  Tandon,  1993, 
Tandon  and  Rosner,  1994) 


Fig.  2  Drag  reduction  associated  with  effective 
permeability  for  a  quasi-spherical  "fractal"  aggregate 
comprised  of  N  primary  spheres  in  the  continuum  regime 
(a  =  Rmax=  1(3/2)  •[Df+2)/Df] }  1/2.Rgyration)  (after 
Tandon  and  Rosner,  1994) 


Fig.  3  Pressure  dependence  ( via  the  Knudsen  number  based 
on  primary  sphere  diameter)  of  the  accessible  surface 
area  of  large  "open"  (Df=1.8)  aggregates;  reaction 

probability,  a,  of  probing  molecule  0.1;  (after  Rosner  and 
Tandon,  1994) 


Fig. 4  "Slof'-type  counterflow  diffusion  flame  (CDF-) 
burner  set-up  for  in  situ  and  extractive  experimental 
studies  of  the  nucleation,  growth,  transport  and 
restructuring  of  aggregates  in  flames  (after  Albagli,  Xing 
and  Rosner,  J994;  see,  also  Gomez  and  Rosner,  1993) 
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Fig.  5  Typical  multiparticle  ^aggregate  thermophoretically 
extracted  from  laminar  CDF  seeded  with  Ti02  precursor 
TOP  vapor.  TEM  image  (a)  compared  to  'touching  sphere’ 
idealization  (b)  (after  Albagli,  Xing  and  Rosner,1994) 
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Fig.  7  Axisymmetric  impinging  jet  CVD-reactor  wit! 

Sm  [0  fte,d  "Postal"  (after  Rosner,  Collins  anc 
illo  1993,  Collins,  1994)  for  systematic  studies  o: 
oxide  film  deposition  from  the  vapor  precursor  TTIP 


r  a  l  ucposHlon  rate  data  (sported  as  an  apparent 
first  order  heterogeneous  rate  constant)  from  TTIP/Cb/N-) 

mainstream  showing  'best-fit'  deposition  rates  (after 
Rosner  Collins  and  Castillo,  1993,  Castillo  and 
Rosner,  1994)  predicted  using  chemical  sublayer  (CSL-) 
theory,  allowing  for  homogeneous  consumption  of  TTIP 
within  thermal  BL. 


2.3.  Multiphase  transport  theory:  nucleation,  growth,  thermophoresis  and 

INERTIA;  AEROSOL  SAMPLING  IMPLICATIONS 

In  this  program  we  have  completed  and  submitted  for  IJHMT  publication  a 
comprehensive  set  of  Seeded  micro-combustor  experiments,  and  ancillary  theoretical 
calculations  on  the  interesting  competition  between  particle  inertia  and  particle 
thermophoresis  for  the  case  of  particle  transport  across  laminar  nonisothermal  gaseous 
boundary  layers  on  surfaces  with  streamwise  curvature  (e.g.,  turbine  blades). 
(Konstandopoulos  and  Rosner,  1994).  For  inertial  impaction  we  demonstrated  that  our 
earlier  idea  of  correlating  impingement  rates  from  compressible  gas  flows  using  an  effective 
Stokes  number  (Israel  and  Rosner,  1983)  can  be  generalized  to  include  the  effects  of 
aerodynamically  interacting  targets,  as  in  a  'cascade'  of  turbine  blades  (Konstandopoulos, 
Labowsky  and  Rosner,  1993). 

We  also  completed  and  published  a  computational  study  of  the  unusual  population 
dynamics  of  coagulating  absorbing-emitting  particles  in  strong  radiation  fields  (Mackowski 
et.al.,  1994).  (For  a  useful  overview  of  our  recent  AFOSR-supported  work  on  these  and 
other  effects  of  energy  transfer  on  suspended  particle  dynamics,  see  Rosner,  et.  ai,  1992). 

The  systematic  effects  of  particle  size/morphology-dependent  wall  deposition  and 
coagulation  on  the  sampling  of  aerosols  have  been  predicted  and  discussed  in  2  papers,  one 
(Rosner  and  Tassopoulos,  1992)  which  appeared  during  this  period,  and  one  now  in 
preparation  (Rosner,  Tandon  and  Konstandopoulos,  1994). 


2.4.  Kinetics  and  morphology  of  cvd-materials  in  multi-phase  environments 

A  small  impinging  jet  (stagnation  flow)  reactor  (Fig. 7)  has  been  developed  and  used 
to  study  the  chemical  vapor  deposition  (CVD-)-rates  of  refractory  layers  on  inductively 
(over-)heated  substrates  (Collins,  Rosner  and  Castillo, 1992, 1993;  Collins,  1994).  These 
measurements,  initiated  with  the  co-sponsorship  of  NASA-Lewis  Labs,  have  been  used  to 
understand  deposition  rates  and  associated  deposit  microstructures  observed  in  highly  non¬ 
isothermal,  often  particle-containing  local  CVD  environments.  Figure  8  shows  (logarithmic 
ordinate)  our  apparent  deposition  rate  probabilities  vs.  reciprocal  surface  temperature  for 
Ti02(s)  obtained  from  TTIP(g).  The  solid  curve  marked  Damhom111  2x  1014  (where  this 
parameter  may  be  regarded  as  a  dimensionless  homogeneous  rate  constant)  shows  our 
predicted  CVD  rate  behavior  including  non-negligible  Soret  transport  (Castillo  and 
Rosner,  1993)  allowing  for  TTIP  decomposition  within  the  boundary  layer.  Agreement  with 
the  experimental  data  of  Collins  (1994)  is  encouraging  indeed.  Regarding  deposit  surface 
morphoplogy  ,  we  completed  a  theoretical  study  of  the  morphological  stability  of  deposits 
growing  by  the  mechanism  of  thermophoresis  (or  Soret  transport  for  vapors)  (chemical 
sublayer  (CSL-)  theory  of  Castillo,  Garcia- Ybarra  and  Rosner,  1992).  On  the  subject  of  the 
'grain'  density  of  vapor  deposits  we  proposed  and  successfully  demonstrated  a  rational 
correlation  which  relates  the  density  of  7  materials  (including  SiC,  BN  and  graphite)  for 
which  data  could  be  found  in  the  literature  to  their  fundamental  properties  (activation  energy 
for  surface  diffusion,  lattice  dimensions,...)  and  deposition  conditions  (  reactant  pressure, 
surface  temperature,  reaction  probability)  via  a  suitable  Damkohler  number  which  we  call 
the  'burial'  parameter  (Kho,  Collins  and  Rosner,  1994). 

In  this  program  we  have  also  summarized  our  research  on  the  high  temperature 
gasification  kinetics  of  solid  boron  by  B2C>3(g),  including  the  implications  of  our  flow  reactor 
data  for  chemical  propulsion  applications  (Gomez,  Rosner,  and  Zvuloni,  1993). 

In  our  OSR-sponsored  Yale  HTCRE  Lab  research  during  this  program,  briefly 
reviewed  here,  we  have  shown  that  new  methods  for  rapidly  measuring  particle  transport 
rates,  combined  with  recent  advances  in  boundary  layer  theory,  provide  useful  means  to 
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identify  and  incorporate  important,  but  often  previously  neglected,  mass  transport 
phenomena  in  many  multiphase  propulsion  engineering  and  materials  engineering 
design/optimization  calculations. 

Despite  the  formidable  complexities  to  be  overcome  in  the  design  and  operation  of 
air-breathing  propulsion  power  plants  utilizing  a  broad  spectrum  energetic  fuels  these 
particular  techniques  and  results  are  indicative  of  the  potentially  useful  simplifications  and 
generalizations  which  have  emerged  from  this  program’s  fundamental  AFOSR-funded 
research  studies  of  combustion-generated  particle  transport  mechanisms.  It  is  hoped  that 
this  Final  Report  and  its  supporting  (cited)  papers  will  facilitate  the  refinement  and/or 
incorporation  of  some  of  the  present  ideas  into  engineering  design  procedures  of  much 
greater  generality  and  reliability.  This  work  has  already  helped  identify  new  directions 
where  research  results  could  have  a  significant  impact  on  engineering  practice  in  both  the 
defense  and  civilian  sectors  of  the  US  economy  (Section  3.4). 


3.  ADMINISTRATIVE  INFORMATION: 

PERSONNEL,  PRESENTATIONS,  APPLICATIONS, 

"COUPLING’  ACTIVITIES 

The  following  sections  summarize  some  pertinent  'non-technical'  facets  of  the 
abovementioned  Yale  HTCRE  Lab/AFOSR  research  program: 

3.1  Personnel 

The  present  results  (Sections  2  and  5)  are  due  to  the  contributions  of  the  individuals 
listed  in  Table  3.1-1,  which  also  indicates  the  role  of  each  researcher  and  the  relevant  time 
interval  of  the  activity.  It  will  be  noted  that,  in  addition  to  the  results  themselves,  this 
program  has  simultaneously  contributed  to  the  research  training  of  a  number  of  students  and 
3  recent  PhDs,  who  will  now  be  in  an  excellent  position  to  make  future  contributions  to 
technologies  oriented  toward  air-breathing  chemical  propulsion,  and  high-tech  materials 
processing. 

Table  3.1-1  Summary  of  Research  Participants a  on  AFOSR  Grant  : 


TRANSPORT  PHENOMENA  AND  INTERFACIAL  KINETICS 
IN  MULTIPHASE  COMBUSTION  SYSTEMS 


Name 

Status^ 

Date(s) 

Principal  Research 

Activity!* 

Albagli,  D. 

PDRA 

4/92-5/94 

particle  prod/char,  in  CDFs 

Cohen,  R.  D. 

VS 

Spring'93 

aggregate  restructuring  theory 

Collins,  J. 

GRA 

'92-'94 

CVD  of  ceramic  coatings 

Gomez,  A. 

Asst.Prof. 

'92 

Ms.  on  particle  transp.  props. (CDFs) 

Kho,T. 

GRA 

'92-'94 

correl.  of  density  of  CVD  coatings 

Konstandopoulos,  A 

GRA, PDRA 

’92,  '93 

combined  inertia  +  thermophoresis 

Labowsky,  M.  J. 

VS 

'91 -'94 

inertial  impaction  and  erosion 

Papadopoulos,  D 

GRA 

'92-'94 

boundary  conditions  at  G/S  interface 

Rosner,  D.E. 

PI 

'91 -'94 

program  direction-dep.  theory/exp 

Silverman,  I. 

PDRA 

'92-'93 

spray  evap/comb.  at  high  pressures 

Tandon,  P. 

GRA 

'92-'94 

transport  phenom.  in  BLs  and  CDFs 

Xing,  Y. 

GRA 

'93-'94 

particle  prod/char,  in  CDFs 

3  PDRA=Post-doctoral  Research  Asst  GRA=  Graduate  Research  Assistant 


PI  =  Principal  Investigator  VS  =  Visiting  Scholar 

b  See  Section  5  for  specific  references  cited  in  text  (Section  2) 


3.2  Cooperation  with  US  Industry 

The  research  summarized  here  was  supported  by  AFOSR  under  Grant  91-0170 
(2/15/91-2/14/94).  The  Yale  HTCRE  Laboratory  has  also  been  the  beneficiary  of  continuing 
smaller  grants  from  U.S.  industrial  corporations,  including  groups  within  GE,  DuPont,  Union 
Carbide  (now  Advanced  Ceramics  Corp.)  and  Shell,  as  well  as  the  feedback  and  advice  of 
principal  scientists/engineers  from  each  of  these  corporations  and  Combustion  Engineering- 
ABB  and  Textron.  We  appreciate  this  level  of  collaboration,  and  expect  that  it  will 
accelerate  inevitable  applications  of  our  results  in  areas  relevant  to  their  technological 
objectives  (see,  also,  Section  3.4,  below). 

3.3  Presentations  and  Research  Training 

Apart  from  the  publications  itemized  in  Section  5  and  our  verbal  presentation  (of 
progress)  at  the  regular  AFOSR  Contractors  Meetings  ,  our  results  have  also  been 
presented  at  some  30  seminars/conferences-including  annual  or  topical  conferences  of  the 
following  professional  organizations: 

Int.  Fine  Particle  Res.  Inst.  (6/92,  6/93) 

AIChE  (11/91,  11/92,  11/93) 

A  AAR  (10/93) 

Electrochemical  Soc.  (CVD  XI  (5/91),  XII;  5/20/93) 

ASME-Engineering  Foundation  (3/91) 

Materials  Research  Society 
Combustion  Inst. 

In  addition,  during  the  period:  2/15/91-2/14/94,  the  PI  presented  seminars  at  the 
following  Universities: 


U  Manchester  5/28/92 
Northwestern  4/22/93 
Brown  (9/22/92) 

U.  Wisconsin 

Trondheim 

Waterloo 


Leeds  5/29/92 
CUNY  10/18/93 
Notre  Dame  (10/27/92) 
U.  Oslo 
U.  Tolouse 


Penn  State  (7/28/92) 
KTH-Stockholm 
U.  Paris-Nord. 

U.  Limoges 
Technion-Haifa 


This  program  involved  completion  of  the  PhD  dissertation  research  of  three  Yale 
graduate  students  (J.  Collins  ,  A.G.  Konstandopoulos  and  M.  Tassopoulos;  cf  Table  3.1-1) 
and  will  form  the  basis  of  the  dissertation  of  P.  Tandon  (to  be  completed  during  '94-'95). 
Also,  T.  Kho  has  received  her  Master's  Degree  based  on  work  supported  in  this  program. 

3.4  Some  Known  Applications  of  Yale-HTCRE  Lab  Research  Results 

It  has  been  particularly  gratifying  to  see  direct  applications  of  some  of  this  generic 
AFOSR-supported  particle  and  vapor  mass  transfer  research  in  more  applications-oriented 
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investigations  reported  in  recent  years.  Indeed,  the  writer  would  appreciate  it  if  further 
examples  known  to  the  reader  can  be  brought  to  his  attention. 

Our  AFOSR  supported  research  on  soot  deposition  rates  from  flowing  laminar  or 
turbulent  combustion  gases  has  been  applied  by  Aerojet  Corp.  (D.  Makel  et.al.,  1990)  to 
develop  a  model  for  application  to  rocket  chambers  and  nozzles  (with  NASA  support). 
Extensions  to  jet  engine  nozzles  are  currently  being  made  by  M.T.  Nys  at  Pratt  and 
Whitney  Engine  Business  in  W.  Palm  Beach  FL 

In  the  area  of  multicomponent  vapor  deposition  in  combustion  systems  applications 
of  our  predictive  methods  (for  "chemically  frozen"  (Rosner  et.al.,  1979)  and  LTCE 
multicomponent  laminar  boundary  layers)  have  been  made  by  British  Coal  Corporation- 
Power  Generation  Branch  (I.  Fantom,  contact)  in  connection  with  their  topping  cycles  which 
run  gas  turbines  on  the  products  of  fluidized  bed  coal  combustors/gasifiers.  Also,  in 
combustion  research  many  groups  (eg.,  Dobbins  et.al.  (Brown  U.),  Faeth  et.al.  (U.  Mich.), 
Katz  et  al.  (J.  Hopkins  U.))  are  now  utilizing  "thermophoretic  sampling"  techniques  to 
exploit  the  size-  and  morphology-insensitive  capture  efficiency  characteristics  that  we  have 
proven  in  our  AFOSR  research  (Section  2.1). 

Our  AFOSR  and  NASA  fundamental  research  on  chemical  element  segregation  in 
the  CVD  of  refractory  ceramics  (eg.,  SiC  and  metal  borides)  (see,  eg.,  Collins  and  Rosner, 
1991,  1992)  is  evidently  of  use  to  AFML  contractors  synthesizing  controlled  stoichiometry 
fibers  for  light  weight/high  strength  composites  (Americom,  Textron). 

For  calculating  suspended  particle  concentrations  along  trajectories  outside  of 
aircraft  (involved  in  atmospheric  sampling),  or  inside  of  CVD  reactors,  A.  S.  Geller  and  D.  J. 
Rader  of  Sandia-Albuquerque  have  adopted  a  method  developed  in  our  earlier  AFOSR 
work  (Fernandez  de  la  Mora,  1981),  and  recently  applied  in  our  own  studies  of  particle 
motion  in  laminar  boundary  layers  with  streamwise  curvature  (Konstandopoulos  and 
Rosnert,  1994). 

Ongoing  work  at  MIT  (Walsh  et.al.  1992),  PSI  (J.J.Helble)  and  Sandia  CRF  (L.L. 
Baxter)  has  incorporated  our  rational  correlation  of  inertial  particle  impaction  (e.g.  a 
cylinder  in  cross-flow)  in  terms  of  our  effective  Stokes  number  (Israel  and  Rosner,  1983, 
and  Konstandopoulos  et.  al.  1993).  Recent  applications  of  our  AFOSR  and  DOE-supported 
research  (on  the  correlation  of  inertial  impaction  by  cylinders  in  crossflow)  have  also  been 
made  by  the  National  Engineering  Laboratory  (NEL)  of  Glasgow  Scotland  (Contact:  Dr.  A. 
Jenkins).  NEL  is  apparently  developing  mass-transfer  prediction  methods  applicable  to 
waste-heat  recovery  systems  in  incinerators,  as  well  as  pulverized  coal-fired  boilers.  These 
applications  are  somewhat  similar  to  those  reported  by  the  Combustion  R&D  group  at  MIT 
and  Penn  State  U,  and  are  now  being  taken  up  by  VTT-Energy /Aerosol  Technology  Group, 
in  Finland. 


Explicit  use  of  our  studies  of  self-regulated  "capture"  of  incident  impacting  particles 
(Rosner  and  Nagaragan,  1987)  is  being  made  in  current  work  on  impact  separators  and 
ceramic  heat  exchangers  for  coal-fired  turbine  systems  in  high  performance  stationary  power 
plants.  Other  potential  applications  arise  in  connection  with  "candle  filters"  used  to  remove 
fines  (sorbent  particles,...)  upstream  of  the  turbines.  A  useful  summary  of  work  in  these 
interrelated  areas  (Solar  Turbines,  Textron  Defense  Systems,  Hague  International,...)  was 
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presented  at  the  Engineering  Foundation  Conference  Inorganic  Transformations  and  Ash 
Deposition  During  Combustion.,  the  proceedings  of  which  appeared  in  1992. 

Clearly,  fruitful  opportunities  for  the  application  of  our  recent  "non-Brownian" 
convective  mass  transfer  research  now  exist  in  many  of  the  programs  currently  supported  by 
the  US  Air  Force,  as  well  as  civilian  sector  R&D. 

4.  CONCLUSIONS 

In  the  OSR-sponsored  Yale  HTCRE  Lab  research  during  the  period:  2/15/91-2/14/94, 
briefly  described  above,  we  have  shown  that  new  methods  for  rapidly  measuring  particle- 
mass  transfer  rates,  combined  with  our  recent  advances  in  mass  transport  theory,  provide 
useful  means  to  identify  and  incorporate  important,  but  previously  neglected,  mass  transport 
phenomena  in  many  chemical  propulsion  engineering  and  materials  engineering 
design/optimization  calculations.  One  important  class  of  examples  involve  our  treatment  of 
aggregated  particle  transport  phenomena  (Section  2.1) 

Despite  formidable  complexities  to  be  overcome  in  the  design  and  operation  of  mobile 
and  stationary  power  plants  utilizing  a  broad  spectrum  of  energetic  fuels  the 
abovementioned  techniques  and  results  (Section  2)  are  indicative  of  the  potentially  useful 
simplifications  and  generalizations  emerging  from  our  present  fundamental  AFOSR-funded 
research  studies  of  combustion-generated  particle  transport  mechanisms  and  interfacial 
reactions  relevant  to  the  synthesis  of  refractory  materials.  It  is  hoped  that  this  Final  Report 
and  its  supporting  papers  (Section  5)  will  facilitate  the  incorporation  of  many  of  the  present 
ideas  into  design  and  test  procedures  of  greater  generality  and  reliability.  This  work  has 
also  helped  identify  new  directions  where  it  is  anticipated  that  research  results  from  this 
AFOSR  program  have  a  significant  impact  on  future  DOD  and  civilian  sector  engineering 
practice. 
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Yale  University,  Department  of  Chemical  Engineering 
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Literature  CVD-rate  data,  and  our  recent  experiments  on  Ti02(s) 
film  growth  from  titanium  (IV)  tetra-isopropoxide  (TTIP)  vapor 
using  a  well-defined  impinging  jet  reactor,  reveal  that  the  onset 
of  vapor  phase  reactions  near  a  hot  deposition  surface  can  lead  to 
sharp  reductions  in  CVD-rate  and  alterations  in  deposit 
microstructure.  These  observations  have  motivated  our 
development  of  a  thin  chemical  sublayer  (CSL-)  theory  for 
predicting  (using  simple  formulae)  the  interplay  of 
heterogeneous  kinetics,  homogeneous  kinetics  and  (Fick-,  Soret-, 
convective-)  transport  phenomena  in  CVD  reactors.  CSL  theory, 
only  briefly  outlined  here,  can  be  used  to  interrelate  different 
CVD  reactors /conditions  and  thereby  guide  the  selection  of 
reactor  configurations/conditions  that  will  lead  to  the 
maximum  CVD-rate  prior  to  the  onset  of  vapor  phase  reactions 
and/or  external  transport  limitations.  Comparisons  with  our 
present  TiC>2(s)  (via  TTIP)  data  are  used  to  suggest  fruitful 
extensions  of  this  work. 


1.  INTRODUCTION 

To  interpret  our  recent  experimental  observations  on  the  CVD  rate  of 
Ti02  films  in  an  impinging  jet  reactor  (1,5)  and  guide  the  design/operation  of 
future  CVD  reactors,  a  rational  asymptotic  theory  has  been  developed  for  the 
onset  of  homogeneous  chemical  reactions  in  the  thermal  boundary  layer  (BL) 
adjacent  to  a  hot  deposition  surface,  with  emphasis  on  their  effect  on 
deposition  rates  and  film  quality  (1,2).  The  analysis  is  tailored  to  systems  in 
which  a)  gas  phase  reactions  are  confined  to  a  thin  chemical  sublayer  (CSL) 
embedded  inside  the  thermal  BL  and  b)  the  homogeneous  reaction(s)  of 
interest  can  be  represented  by  a  single  high  activation  energy  (Ehom/(RTe)»  1) 
Arrhenius-type  rate  expression.  In  this  limit  one  finds  (2): 
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ScSL  ^(EhomV1.,^ 

8h  lRTe/  0W-1  [1] 

where  6W  (=  Tw/Te)  is  the  temperature  ratio  across  the  thermal  BL  of  slope 

thickness  8h-  Because  of  the  large  molecular  weight  disparities  and  high 
temperature  gradients  prevailing  in  many  CVD  systems  (3,4),  we  allow  for  the 
Soret  reduction  of  dilute  reagent  Fick  transport  to  the  hot  surface.  The 
resulting  closed-form  expressions  provide  CVD-rate  predictions  and  rational 
quantitative  criteria  for  "vapor  phase  ignition"  (VPI)  in  terms  of  CVD  system 
parameters  and  the  "known"  chemical  kinetic  parameters  characterizing  the 
vapor  reactants.  Conversely,  armed  with  such  a  theory,  one  can  use 
experimentally  observed  VPI  conditions  to  infer  the  effective  homogeneous 
kinetic  parameters  for  the  system  in  question— information  often  not 
independently  available.  In  effect,  we  are  recommending  use  of  the  gaseous 
boundary  layer  itself  as  a  "flow  reactor",  with  the  CVD  surface  acting  as  a 
detector  of  the  "remaining  reactant."  Our  general  results  can  be  easily  applied 
to  specific  ceramic  film  systems  of  current  interest,  with  our  present  emphasis 
being  the  CVD  of  Ti02(s)  (from  TiCU  or  TTIP+  O2) 

The  overall  trends  predicted  by  CSL  theory  are  displayed  in  Fig.  1  over  a 

broad  range  of  heterogeneous  Damkohler  numbers  Damhet  0s  Ahet  5/D)- 
values  and  Te/Tw-values,  for  particular  dimensionless  activation  energies 

Ehom/(RTe)  and  Ehet/(RTe),  homogeneous  Damkohler  number  Darnhom  (s 

Ahom  52/D)  and  Soret  factor  <xT  (2,3).  Here  Ahet  and  Ahom  are,  respectively,  the 
heterogeneous  and  homogeneous  rate  constant  pre-exponential  factors.  In  the 
CVD  reactor  described  below(l,5)  using  TiCU(g)  as  the  Ti-carrier  we  were 
unable  to  study  such  VPI  trends  at  convenient  surface  temperatures, 
however,  the  use  of  TTIP  provided  what  appears  to  be  a  clear  cut  example  of 
this  phenomenon  (see  Fig.  3  below).  From  available  estimates  of  the  overall 
homogeneous  chemical  kinetics  of  this  system,  we  infer  Ehom/  (RTe)  ®  22  and 
0W  *  3,  which  implies  (Eq.[l])  that  ScsiVSh  -  1/5.  This  thickness  ratio  may  be 
small  enough  to  apply  chemical  sublayer  theory  as  a  useful  first 
approximation,  but  because  of  its  magnitude  we  have  also  initiated  a  finite- 
element  numerical  attack  on  this  same  example  (8),  as  well  as  needed 
generalizations  of  CSL-theory.  Incidentally,  an  interesting  previously 
observed  example  of  VPI  using  TiCLi(g)  (6)  could  not  be  analyzed  using  our 
present  approach  due  to  uncertain  (undocumented)  transport  conditions  . 

2.  EXPERIMENTAL 

We  are  now  investigating  "vapor  phase  ignition"  (VPI)  under  well- 
defined  transport  conditions  using  an  axisymmetric  impinging  laminar  jet 
"pedestal”  CVD  reactor  (Fig.  2).  Upon  entering  the  reactor,  cold  (300-400  K) 


At 


reagent  vapor  (TiCU  or  TTIP)  and  carrier  gas  (Ar  or  N2)  flow  through  a  short 
mixing  chamber  and  a  converging  nozzle.  The  laminar  jet  emerging  from  the 
nozzle  impinges  on  a  polished  quartz  substrate  (dsubstrate  =  djet  =  12.7mm; 
located  roughly  one  jet  diameter  downstream)  heated  from  below  by  an  RF- 
heated  graphite  susceptor  .  CVD-rate  measurements  are  by  in  situ 
interferometry  (at  633  nm)  and/or  ex  situ  weight  gain.  Deposit 
microstructures  are  studied  using  scanning  electron  microscopy  (SEM). 

Preliminary  experiments  were  for  titania  deposition  from  TiCl4/C>2/Ar 
and  TiCl4/N20/Ar  at  0.1  MPa  total  pressure(l).  To  explore  VPI  phenomena  at 
more  accessible  substrate  temperatures,  we  shifted  to  titanium(IV)-iso- 
propoxide  (Ti(OC3H7)4)  in  N2  with  a  small  amount  of  added  oxygen  to  "bum 
away"  any  co-deposited  carbon.  The  liquid  reagent  source  is  a  constant 
temperature  bubbler.  All  reagent  lines  are  teflon  or  stainless  steel,  heated 
above  the  reagent  vapor  dew  point.  To  avoid  water  contamination  the  reactor 
is  pumped  down  to  10  Pa  for  several  hours  before  each  run  and  only  ultra 
high  purity  Ar  and  N2  gases  are  used  ([H2O]  <  0.5  ppm). 

Our  current  range  of  accessible  operating  conditions  (with  well-defined 
fluid  flow)  are:  substrate  temperature  Tw  up  to  1600  K,  total  pressure  p 
between  0.01  and  0.1  MPa,  and  impinging  jet  Reynolds  number,  Rejet , 
between  roughly  200  and  800  (based  on  nozzle  diameter,  djet).  Under  these 
conditions  natural  convection  is  not  expected  in  the  jet  impingement  region 

since  Rej^/Gij/4  >>1,  (Grh  is  the  heat  transfer  Grashof  number  for  the 
impinging  jet  boundary  layer  (see,  e.g.,( 9))).  Moreover,  we  estimate  negligible 
entrainment  of  recirculating  reaction  products  in  the  impingement  region. 

3.  RESULTS 

As  shown  in  Fig.  3,  in  the  CVD  of  TiC>2  from  TTIP  above  Tw  -  1150K  we 
have  observed  a  deposition  rate  fall-off  to  about  l/10th  of  the  diffusion- 
controlled  plateau  value.  In  Fig  3,  keff  is  an  effective  first  order  rate  constant 
defined  as  the  absolute  molar  deposition  flux  (-jw"/M)rrip  divided  by  a 
reference  concentration  [Pmp,e/(RTW)].  The  apparent  activation  energy  of  the 
heterogeneous  kinetics  controlled  region  below  750  K  is  about  132  KJ/mole. 
Before  reaching  VPI,  our  data  go  through  a  broad  diffusion-controlled 
plateau.  The  actual  rate  fall-off  region  above  ca.  1150K  shows  a  slope  of 
roughly  -70  KJ/mole.  The  break-point  temperature  may  be  influenced  by 
reagent  concentration,  total  pressure  and  flow  rate,  and  we  are  currently 
investigating  these  possibilities.  Preliminary  analysis  of  SEM  micrographs  of 
Ti02  deposits  from  TTIP  shows  that  relatively  smooth,  dense  polycrystalline 
films  are  produced  under  heterogeneous  kinetics-controlled  conditions  (Tw  < 

750K).  Under  diffusion-controlled  conditions  we  observe  massive  (ca.  40pm) 
and  poorly  adherent  dendritic  structures.  Interestingly,  films  produced  well 
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above  the  onset  of  VPI  show  higher  density  and  improved  substrate 
adherence. 

4.  DISCUSSION 

Using  CSL-theory(2)  we  can  readily  predict  the  wall  concentration  and 
deposition  flux  of  reagent,  allowing  for  depletion  due  to  CSL  reactions.  In 
cases  where  the  homogeneous  reactions  generate  non-depositing  products, 
this  'surviving'  reagent  deposition  flux  should  equal  the  experimentally 
observed  deposition  rate.  The  closed  form  asymptotic  theory  yields  all 
expected  transitions  in  deposition  regimes  for  CVD  systems  in  which  Ehom 
>Ehet  (cf.  Fig.  1).  At  sufficiently  low  surface  temperatures  the  homogeneous 
reaction  is  suppressed  and  the  CVD-rate  is  limited  by  the  intrinsic 
heterogeneous  chemical  kinetics.  However,  as  Tw  increases,  gas  phase 
transport  limitations  set  in  (see  below)  and  TTIP  vapor  reactions  become 
increasingly  important,  ultimately  causing  the  deposition  flux  of  virgin 
reagent  to  decrease  (Figs. 1,3).  In  the  present  (TTIP/O2)  system,  this  occurred 
after  the  broad  external  diffusion-controlled  plateau  (750  <  Tw  <  1150K) 
revealed  in  Fig.3. 

To  within  our  experimental  errors  (±10%)  the  'plateau'-value  of  keff  is 

seen  to  correspond  to  our  estimate  of  6h,eff  based  on  earlier  impinging 
laminar  jet  heat  transfer  measurements,  provided  a  significant  correction  is 
made  for  the  Soret-reduction(3,4)  in  diffusion-controlled  TTIP  mass  transfer 

rate  ( ca .  Vs  at  VPI,  based  on  an  estimated  ax- value  of  1.6)  and  systematic 
corrections  for  non-constant  properties  (ca.  15%)  and  nonunity  Lewis  number 

(ca.  33%,  since  DTTlP-N2/(k/(pcp))  -  0.26).  No  appreciable  Stefan  flow 
("suction")  correction  was  necessary  (see,  eg.,(10))  since,  in  these  experiments 
the  mass  fraction  of  TTIP  in  the  feed  was  only  0.5%,  corresponding  to  a 
titanium  element  mass  fraction(9)  of  only  0.08%.  Indeed,  the  present  version 
of  CSL  theory(2)  also  neglects  Stefan  flow  phenomena,  as  is  appropriate  in  the 
'dilute'  reagent  limit. 

For  a  first  order  homogeneous  reaction  CSL-theory  predicts  that  VPI 
should  occur  approximately  where: 

0w,VPI  =  ( Ehom/(RTe))*(  ln[  Damhom*(Ehom/(RTe))-1]}-1  [2] 

when  VPI  occurs  along  the  diffusion-controlled  branch  (Figs.1,3).  The 
experimentally  observed  value,  8W,VPI  *1150/380  =3.0,  is  in  approximate 
accord  with  published  values(ll)  of  Ehom  (70.5  KJ/mole)  and  Ahom(4xl05  s'1), 
since,  under  our  experimental  CVD  conditions  we  expect  Damhom“O(104)* 
While  CSL  theory(2)  also  predicts  that  the  ultimate  slope  of  the  fall-off  region 
should  correspond  to  an  apparent  activation  energy  of  -(Ehom  -  Ehet)/  the 
observed  slope  appears  closer  to  -Ehom  •  This  may  indicate  that  at  such  high 
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surface  temperatures  Ehet  has  decreased  to  values  small  compared  to  Ehom- 
Extended  CVD-rate  data  and  future  generalizations  of  CSL-theory  should 
clarify  this  situation,  and  provide  further  valuable  tests  of  CSL  theory. 

Returning  to  the  onset  of  appreciable  external  diffusion  limitations,  for 
a  first  order  heterogeneous  reaction  this  should  occur  at  the  Tw-value  for 
which  Damhet*exp["  Ehet/(RTW)  ]  -1.  This  point,  also  shown  on  Fig.3,  is 
consistent  with  our  own  lower  temperature  (Ahet/  Eher)date.  In  considering 
the  two  "end-points"  shown  on  Fig.  3,  we  remark  that  in  many  situations 
merging  them,  if  acheivable,  would  be  an  optimum  operating  point.  This 
would  maximize  chemically-controlled  CVD-rates  while  just  avoiding 
increased  sensitivity  to  position  and  gas-dynamic  conditions  (leading,  also,  to 
surface  roughness),  as  well  as  avoiding  deposit  defects  associated  with  the 
capture/incorporation  of  embryonic  microparticles. 

5.  CONCLUSIONS  and  RECOMMENDATIONS 

Using  Ti02  CVD  rate  data  (1,5,  Fig.  3),  we  have  demonstrated  that  an 
asymptotic  chemically  reacting  sublayer  (CSL-)  theory(2)  can  provide  useful 
estimates  of  unreacted  reagent  concentrations  and  wall  fluxes  at  CVD 
surfaces.  This  opens  the  way  to  rationally  interrelate  different  CVD  reactors 
environments  via  the  intermediary  of  their  global  reaction  kinetic 
parameters  (2).  Simple  CSL-theory  is  applicable  and  convenient  for  reactor 
design/optimization  whenever  homogeneous  reactions  can  be  represented  by 
a  global  reaction  with  Eh0m/(RTe)»  1  and  the  temperature  ratio  across  the 
boundary  layer  is  also  sufficiently  large  (c/.  Eq.  [1]).  This  approach  could  also  be 
exploited  in  computational  methods  when  high  activation  energies  and  large 
local  gradients  preclude  accurate  numerical  treatment  of  the  actual  network 
of  homogeneous  chemical  reactions.  Thus,  the  structure  of  such  "chemical 
sublayers"  can  be  computed  (using  "rescaled”  coordinates)  in  greater  chemical 
detail  by  exploiting  geometrical  and  property  value  simplifications  valid 
because  of  their  thinness.  Indeed,  even  when  conventional  boundary  layer 
theory  breaks  down  (due  to  inadequate  (Re-Pr)1/2;  see,  eg,  (9)),  boundary  layer 
approximations  may  remain  valid  and  useful  for  such  chemical  sublayers  — 
as  will  be  demonstrated  our  follow-on  research. 

Our  current  investigations  include  obtaining  direct  experimental 
corroboration  of  vapor  phase  reaction  onset,  generalizing  CSL-theory,  and 
developing  a  self-consistent  theory  which  explicitly  accounts  for  the 
production,  growth  and  transport  of  clusters  within  CVD-boundary  layers. 
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Fig.  1:  CSL-theory  predicted  Arrhenius 
behavior  for  CVD-rates  in  the  presence 
(at  high  surface  temperatures)  of  vapor 
phase  reactions  yielding  non¬ 
depositing  products  (2);  Dependence  of 
nondimensional  CVD-rate  (log  scale) 
on  (reciprocal)  surface  temperature 
and  heterogeneous  reaction 

Damkohler  number  AhetS/D,  where 
Ahet  is  the  pre-exponential  factor  of 
the  first-order  heterogeneous  rate 
constant. 


Fig.  2:  Axisymmetric  impinging  jet 
CVD  reactor  with  inductively-heated 
"pedestal";  after  Collins  (1991,  1993) 
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Fig.  3:  TiC>2  Deposition  rate  data 
from  TTIP/O2/N2  with  calculated 
diffusion-controlled  limits  showing 
effect  of  Soret  diffusion 


REPORT  DOCUMENTATION  PAGE 


form  Approved 
OMf  NO.  070*4  tit 


FvWk  nwnw  >■*»»  •»  P»  tfWirw  ¥  n**""***"  ■»  «wy»*  *»  _»-  J,l'^t-ni|,n  ¥  iKtorwuu ««  t*«P  (»«**"■  r~3S^iN»  t>r#w  wgyt  »>  »"i  g,**'.f??*^l.fl  °** 


1.  agency  usi  omit  (if*  *•*)  p  ,tPO"T  P*T»  ^  Y#5frwt[^ 4*W ScTj? (+)*7 h  7^3 

I  — -  - -  5T  rUNOMO  NWMHM 

«.  mu  amd  surrrru 

INERTIAL  DEPOSITION  OF  PARTICLES  FROM  POTENTIAL  PE  -  61102F 
FLOWS  PAST  CYLINDER  ARRAYS  PR  -  2308 

_ _ _ — -  SA  -  BS 

A-AUTHOAtf)  G  -  2)-fll70(ArD5R) 

Atiianasios  G.  Konstandopoulos,  M.C..AEL  J.  Labowsky  and  Dan.el  E.  Rosner 


7  rUVOKMMG  ORGANUATION  NAME(S)  AND  ADOAESSUS) 

HIGH  TEMPERATURE  CHEMICAL  REACTION 
ENGINEERING  LABORATORY 

YALE  UNIVERSITY 
BOX  200/80  YALE  STATION 
NEW  HAVEN,  CONNECTICUT  06520  U.S.A. 

».  jrOKJOAJNC/MONriOWNG  AGENCY  NAMt(S)  AND  ADOMSSttS) 

AFOSR/NA 
Building  A 1 0 

Bolling  AFB  DC  20332-6448 


■t  PEMORMING  ORGANIZATION 
R I PORT  NUNMER 


10.  SPONSORING / MONITORING 
AGENCY  REPORT  NUMBER 


11.  SUPPLEMENTARY  NOTH 


12».  DISTRIBUTION  /  AVAILABILITY  IT  AT  I  MINT 


Approved  for  public  release;  distribution  is 
uni imited 


13.  ABSTRACT  (NUjumum  200  orordt) 


12b.  DISTRIBUTION  COOl 


Abstract — Wc  exploit  the  computationally  efficient  method  of  images  (MOI)  to  create  steady 
potential  flows  past  finite  and  infinite  circular  cylinder  arrays,  in  order  to  study  inertial  impaction  of 
particles  on  them.  Although  such  flows  are  often  encountered  in  practice  (e.g.  high  Reynolds  number 
“dusty  gas'*  flows  past  heat  exchanger  tube  banks  and  filter  screens),  little  is  yet  known  about  the 
proximity  effects  of  adjacent  collectors  on  the  capture  efficiency  of  any  particular  collector  in  the 
array.  Our  results  for  linear,  symmetric  arrays  in  steady  cross-flow,  containing  an  odd  number  of 
collectors  (from  three  to  infinity),  over  a  wide  range  or  array  spacings,  support  the  conjecture  that  an 
appropriately  defined  erectil  e  Stokes  number  based  on  the  stagnation  region  fluid  deceleration  rate 
(computed  and  reported  here  as  a  function  of  the  number  of  collectors  and  their  spacing)  can 
adequately  correlate,  for  practical  purposes,  the  capture  efficiency  or  a  representative  collector  in  the 
array.  This  scheme,  combined  with  well-known  correlations  for  an  isolated  collector,  simplifies 
significantly  the  task  of  estimating  inertial  impaction  rates. 


14.  SUBJECT  TERMS 


inertial  irnfact/an; ,  c/lbders  in  crossflow 


is.  NllMttt  or  PAGES 

1*3 


14.  PRICE  coo* 


ii'  »OI„T,  CU-iWlCAI^  pr  Hfij.lt.  CLMWIWOT  |l»  y  £V^7',,ICA"0"  r  "°W 

Unclass  J.f  led  Unclassified  I  Unclassified  |  UL 


SUndRid  Fwm  298  (890104  DrifTj 

>ika»<  h  ami  hr  w-u 

HMI 


NSN  7540-01-260-5500 


J.  Aeroutf  Set. 
Printed  tn  Great 


2i(4)A7l-48Z(/}93) 

it  Britain  / 


0021-8502  -93  $6.00  +  0  00 
f  •  1993  Pergamon  Pms  Ltd 
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Abstract — We  exploit  the  computationally  efficient  method  of  images  (MOl)  to  create  steady 
potential  flows  past  finite  and  infinite  circular  cylinder  arrays,  in  order  to  study  inertial  impaction  of 
particles  on  them.  Although  such  flows  are  often  encountered  in  practice  (e.g.  high  Reynolds  number 
“dusty  gas"  flows  past  heat  exchanger  tube  banks  and  filter  screens),  little  is  yet  known  about  the 
proximity  effects  of  adjacent  collectors  on  the  capture  efficiency  of  any  particular  collector  in  the 
array.  Our  results  for  linear,  symmetric  arrays  in  steady  cross-flow,  containing  an  odd  number  of 
collectors  (from  three  to  infinity),  over  a  wide  range  of  array  spacings,  support  the  conjecture  that  an 
appropriately  defined  effective  Stokes  number  based  on  the  stagnation  region  fluid  deceleration  rate 
(computed  and  reported  here  as  a  function  of  the  number  of  collectors  and  their  spacing)  can 
adequately  correlate,  for  practical  purposes,  the  capture  efficiency  of  a  representative  collector  in  the 
array.  This  scheme,  combined  with  well-known  correlations  for  an  isolated  collector,  simplifies 
significantly  the  task  of  estimating  inertia!  impaction  rates. 


NOMENCLATURE 


«. 

n‘ 

R 

S 

Stk 


u 

“p 

U, 

x 

y 


particle  diameter 

x-direction  unit  vector 

outward  unit  normal  of  the  ilh  cylinder 

cylinder  radius 

array  spacing 

Stokes  number 

time 

fluid  velocity  vector 
particle  velocity  vector 

fluid  velocity  upstream  of  the  collector  (at  —  oo) 

coordinate 

coordinate 


Greek  letters 

p  doublet  strength 
p,  dynamic  viscosity  *T 
pp  particle  mass  concentration  (density) 
i  particle  relaxation  time 
<1>  fluid  velocity  potential  # 

i  impact 

sp  stagnation  point 
w  at  the  wall 
oc  at  infinity 


Subscripts 


1.  INTRODUCTION 

Particle  inertial  behavior  results  from  momentum  non-equilibrium  between  the  suspended 
particles  and  the  host  flow  and  is  characterized  by  the  dimensionless  Stokes  number,  Stk 
(ratio  of  particle  momentum  relaxation  time,  t  to  a  characteristic  flow  time-scale,  Trtow): 

Stk  =  —  (1) 

•  Based  on  Paper  5.E6  presented  at  the  1990  AAAR  Annual  Meeting,  18-22  June.  Philadelphia,  PA,  U.S.A. 

*  Present  address:  45  Tselcpi  Str.,  Thessaloniki.  GR-543  52,  Greece. 
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where  (cf.  Nomenclature): 


xJf’C 

(2) 

18/t, 

R 

low  >1 

00 

(3) 

with  pt  being  the  intrinsic  particle  density,  dp  its  diameter,  pt  the  dynamic  viscosity  of  the 
carrier  gas  and  C  the  Stokes-Cunningham  slip  correction  factor  accounting  for  deviations 
from  continuum  drag  behavior  (see, e.g.  friedlander,  1977). 

The  important  role  of  inertia  in  transporting  large-sized  aerosol  particles  (more  specific¬ 
ally  particles  with  Stk-values  larger  than  a  critical,  flow-field  dependent  value,  Stkctil)  to 
collector  surfaces  leading  to  direct  impaction  has  been  recognized  long  ago  (see  Appendix) 
and  the  study  of  Stk  >  StkcrlI  (supercritical)  inertial  effects  is  a  mature  subject  in  aerosol 
science  (e.g.  Fuchs,  1964;  Davies,  1966;  Friedlander,  1977).  Particle  deposition  rates  are 
usually  expressed  in  terms  of  the  dimensionless  capture  fraction,  tjc,p  e[0,  1],  defined  as  the 
ratio  of  the  actual  flow  area  (at  upstream  infinity),  where  depositing  particle  trajectories 
originate,  over  the  projected  collector  area.  Too  numerous  calculations  of  inertial  capture 
efficiencies  exist  in  the  literature  for  well-known  analytically  or  computationally  generated 
fluid  flow  fields  around  isolated  collectors  to  be  explicitly  mentioned  here.  A  particular 
geometry,  that  of  an  isolated  cylinder  in  cross-flow,  has  received  considerable  attention  due 
to  its  relevance  to  several  practical  applications  (e.g7filtration  by  fibrous  filters,  fouling  of 
heat  exchanger  tube  banks,  etc.). 

Inertial  impaction  of  particles  on  an  isolated  cylinder  target,  especially  at  high  Reynolds 
numbers,  Re  (when  potential  flow  provides  a  good  approximation  of  the  flow  field  away 
from  the  surface  of  the  collector)  is  now  rather  well  understood.  Classic  studies  in  the 
aerosol  literature  (for  detailed  references  see  Fuchs,  1964;  Friedlander,  1977)  have  been 
extended  to  cover  large  parameter  spaces  and  useful  correlations  for  important  quantities  of 
interest  have  become  available  for  routine  engineering  predictions.  In  particular,  Israel  and 
Rosner  ( 1 983)  devised  a  correlation  scheme  for  the  inertial  capture  efficiency,  i/c.p  of  isolated 
cylinders  and  spheres  in  steady  potential  flow,  introducing  the  notion  of  an  “effective  Stokes 
number”  Stkcff.  Non-linear  drag  was  accounted  for,1  following  a  suggestion  of  May  and 
Clifford  (1967)  to  use  the  actual  particle  stopping  time  in  the  definition  of  Stkefr.  Differences 
in  collector  geometry  (sphere  vs  cylinder)  were  accommodated  employing  as  the  character¬ 
istic  flow  time  the  respective  flow  deceleration  rale  at  the  stagnation  point  of  each  collector. 
It  was  then  possible  to*collapse  the  characteristic  S-shaped,  inertial  capture  efficiency  curves 
for  the  cylinder  and  the  sphere,  in  terms  of  a  single,  “universal”  curve.  The  scheme  was  later 
applied  by  Wang  (1986)  in  his  study  of  inertial  impaction  and  rebound  of  particles  from 
cylinders  and  by  Wcsscl  and  Righi  (1989)  who  provided  extensive  correlations,  not  only  for 
the  inertial  capture  efficiency  of  a  cylinder  in  potential  flow,  but  also  for  impact  velocities, 
particle  concentrations  at  the  cylinder  surface,  and  angular  distributions  of  quantities  of 
interest.  A  correlation  of  numerical  results  for  ije,p  for  inertial  impaction  from  viscous  flows 
(30  <  Re  <40,000)  past  a  single  cylinder  is  given  in  Ilias  and  Douglas  (1989). 

Inertial  impaction  on  collector  assemblies  is  less  well-studied  than  the  isolated  collector 
cases.  Except  for  some  special  cases,  little  is  yet  known  about  the  “proximity”  effects  of 
adjacent  collectors  on  the  inertial  deposition  rates  especially  from  high  Re  flows  on  any 
particular  cylinder  in  an  array.  This  work  was  motivated  by  the  fact  that  such  situations 
occur  often  in  practice,  e.g.  in  the  entrance  region  to  a  heat  exchanger  tube  bank  or  a  fibrous 
filter  screen,  as.well  as  in  laboratory  experiments  simulating  such  arrangements.  Neighbor¬ 
ing  collector  effects  have  been  studied  more  extensively  for  low  Re  flows  of  immediate 
relevance  to  filtration  processes,  either  in  “unit  cell”  models  (se^  e.g.,Tien,  1989)  or  for  an 
infinite  array  of  fibers  in  cross-flow  (e.g./Tsiang  et  a!.,  1982;  McLaughlin  et  al.,  1986;  Ingham 


*  Low  transonic  How  compressibility  effects  could  also  be  taken  into  account  using  the  Stk,„  technique  (Israel 
and  Rosner,  1983). 
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et  at.,  1989).  Choudhary  and  Gentry  (1977)  apparently  were  the  first  to  study  inertial 
impaction  from  steady  potential  flow  on  a  row  of  two  cylinders,  using  the  method  of  images 
(MOl)  to  construct  the  inviscid,  incompressible  flow  field.  Ingham  el  at.  (1989)  employed  the 
Boundary  Element  Method  (BEM)  to  compute  the  flow  field,  and  focused  on  the  infinite 
array  problem.  These  authors  have  also  considered,  in  addition  to  circular  cylinders, 
elliptical  cylinders.  Particle  trajectories  were  traced  in  these  flow  fields  and  fjcip  of  a  repres¬ 
entative  collector  geometry  was  evaluated  as  a  function  of  Stk.  Results  in  both  studies  were 
given  only  for  a  particular  cylinder  center-to-center  spacing  of  5.35  radii.  Other  more 
elaborate  computations  for  a  cylinder  in  a  “deep”  tube  bank,  using  variants  or  the  krc 
model  of  turbulence  for  the  flow  field  and  incorporating  particle  tracking  modules,  have 
also  appeared  recently  (Schuh  et  at.,  1989;  Fan  el  at.,  1991)  in  the  interest  of  predicting  tube 
erosion.  Both  studies,  however,  presented  only  the  results  of  representative  numerical 
calculations  but  no  »/« P  -(Stk, . . curves. 

In  the  present  work  we  study  inertial  impaction  on  cylinder  arrays  with  the  following 
objectives; 

(1)  To  investigate  how  the  number  of  collectors  and  their  spacing  alter  the  collection 
characteristics  of  a  single  (“test")  cylinder  in  an  array. 

(2)  To  find  appropriate  ways  to  correlate  such  data  and  provide  approximate,  yet 
rational  methodologies  for  practical  situations  that  would  circumvent  the  need  for  further 
detailed  computations  of  the  type  that  we  have  performed. 

We  do  not  address  here  the  possibility  of  particle  rebound,  so  that  in  what  follows 
impaction  and  capture  are  treated  as  synonymous.  The  structure  of  the  present  papier  is  as 
follows:  In  section  2  we  outline  the  numerical  calculation  of  the  flow  field.  In  section  3 
we  provide  the  results  of  our  particle  impaction  computations  on  a  representative  cylinder 
in  finite  and  infinite  collector  arrays.  We  discuss  our  correlation  scheme  for  tjc ,p  in 
section  4  and  apply  it  to  the  data  of  section  3,  while  section  5  summarizes  the  conclusions 
drawn  from  this  work.  An  Appendix  clarifies  some  issues  concerning  the  critical  Stokes 
number  for  particle  impaction  on  bluff  bodies. 


2.  FLOW  FIELD  COMPUTATION 

The  inviscid  flow  field,  u  =  -V<P  past  an  assembly  of  circular  cylinders  is  determined  in 
terms  of  the  velocity  potential  <t>,  which  satisfies  Laplace’s  equation: 

VJ<D  =  0  (4) 

subject  to  the  boundary  conditions  of  zero  normal  velocity  on  all  cylinder  walls: 

V<J>  n'  =  0  i  =  l,2,. .  ,NC„  (5) 

and  a  uniform  flow  at  infinity  (,x  =  -oo): 


Although.in  principle, these  equations  can  be  solved  with  a  variety  of  methods,  some  of 
the  most  efficient  are  those  exploiting  the  superposition  of  fundamental  solutions  of 
Laplace’s  equation  in  such  a  way  as  to  satisfy  the  boundary  conditions  (Fletcher,  1988). 
Superposition  (image)  methods  for  the  solution  of  Laplace’s  equation  are  well  documented 
in  classic  fluid  mechanics  texts  (Milne-Thomson,  1955;  Lamb,  1932).  In  what  follows  we 
therefore  provide  only  an  “algorithmic,”  sketch  of  the  computational  implementation.  While 
our  MOI  code  can  handle  arbitrary  cylinder  arrangements  of  various  radii,  here  we 
concentrate  on  the  practically  important  geometry  of  a  straight  array  of  identical  cylinders 
in  cross  flow.  For  convenience  we  have  used  arrays  containing  an  odd  number  of  cylinders 
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with  the  central  cylinder,  located  at  the  origin  serving  as  our  “test’  collector.  The  arrays  are 
characterized  by  the  number  of  cylinders  Nc„  and  their  dimensionless  center-to-center 
spacing.  S  expressed  in  cylinder  radii  units.  Potential  now  around  an  isolated  cylinder  is 
equivalent  to  the  superposition  of  a  doublet  (dipole)  singularity  (fundamental  solution  of 
Laplace’s  equation  in  2-D)  at  the  axis  of  the  cylinder  and  a  uniform  stream  at  infinity. 

_<J,  =  _((I)<1  +  <l»s)=l/G,KJ^I^I  +  i;«x.  (8) 

Picking  the  velocity  or  the  uniform  How  al  upstream  infinity,  Um,  directed  along  the 
positive  x-axis,  as  our  velocity  scale  and  the  cylinder  radius  R  as  the  length  scale  the 
calculation  of  the  dimensionless  flow  field  proceeds  as  follows: 

(!)  Wc  first  place  a  doublet  at  the  center  of  each  cylinder  in  ihe  assembly  of  unit  strength 
(i.e.  normalized  with  l/*.R2)  aligned  with  the  negative  x-axis.  We  term  these  the  first 

qeneration  doublets.  ....  .  •  r  „n 

(2)  Using  the  so-called  image  relations  we  compute,  for  each  cylinder,  the  images  or  all 

doublets  external  to  the  cylinder.  The  image  of  a  doublet  of  dimensionless  strength 
u  positioned  at  a  dimensionless  distance  (  from  a  cylinder  center  has  a  dimensionless 
strength  «//2  and  is  located  at  the  inverse  point  of  the  line  connecting  the  cylinder  center 
and  the  doublet,  i.e.  at  a  dimensionless  distance  or  1//  from  the  cylinder  center  with  an 
orientation  antiparallel  to  that  of  the  doublet.  These  doublets  form  the  second  generation. 

(3)  We  repeat  the  image  generation  process  of  item  2  for  the  second  generation  doube  s 
creating  a  third  generation  of  doublets  and  so  on,  checking  at  each  iteration  if  the  BC  or 
zero  velocity  normal  to  all  cylinder  surfaces*  (impenetrability)  is  satisfied  to  within  a  pre¬ 
scribed  tolerance  (typically  10-*).  At  this  point  we  have  a  distribution  of  doublets  which 
when  taken  together  with  the  uniform  flow  al  upstream  infinity  satisfy  Laplace  s  equation 
for  the  velocity  potential  and  the  associated  BCs.  For  the  irrotational  flow  considered  here 
this  distribution  of  doublets  suffices  to  determine  uniquely  the  flow  field. 

The  MOI  is  computationally  very  simple  and  typical  programs  have  relatively  few  lines  of 
code.  The  above  mentioned  loop  is  used  repeatedly  to  calculate  new  generations  of  doublets 
by  applying  the  image  relations  to  the  previous  generations  of  doublets.  One  major 
advantage  or  the  MOI  over  other  techniques  is  the  ability  to  calculate  the  velocity  field  at 
any  point  in  the  flow  and  not  just  at  selected  grid  or  node  points.  In  this  way  interpolation 
dThe  flow  field  during  the  calculation  or  particle  trajectories  is  avoided  and  thus  a  source  ol 

numerical  error  is  eliminated.  , 

The  number  of  doublets  generated  with  the  MOI,  Afd,  depends  on  N?,  and  the  number  or 

doublet  generations  performed.  After  N  generations  we  find  that  Nd  is  given  by. 


1-1 


(9) 
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As  may  be  seen  from  Fig.  1,  the  number  or  doublets  grows  exponentially  with  the  number  of 
generations  and  the  number  of  cylinders  in  the  array.  The  number  of  generations  required 
to  obtain  an  accurate  result  depends  on  the  array  geometry.  For  close  spacings  (S<2.5) 
and/or  many  collectors,  the  CPU  memory  (and  time)  requirements  of  the  MOI  increase 
substantially.  For  example,  five  generations  for  an  array  containing  15  cylinders  are 
sufficient  to  approach  the  memory  limit  of  the  workstation  used  in  these  calculations.  We 
have  therefore  explored  several  series  acceleration  techniques  using  the  results  of  intermedi¬ 
ate  steps  to  improve  convergence.  Or  these  techniques,  Shank’s  transformation  (see,e.g.. 
Van  Dyke  1974)  was  found  to  be  well  suited  for  our  calculations.  For  example,  convergence 
to  the  fifth  decimal  digit  could  be  obtained  for  cylinder  assemblies  containing  up  to  five 
cylinders  after  five  doublet  generations  and  repeated  application  or  Shank’s  transformation. 

•We  found  that  monitoring  the  convergence  of  the  centra!  cylinder's  stagnation  point  velocity 
was  a  conservative  criterion  for  satisfaction  of  the  BCs  on  all  cylmder  surfaces.  I»  1k«  ^  cyll 
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Fig.  1.  Dependence  of  number  of  doublets  on  number  of  generations  and  number  of  collectors  in 

the  array( 


This  is  equivalent  to  the  results  obtained  after  10  generations  without  acceleration.  Based 
on  these  results  Shank’s  transformation  was  permanently  incorporated  in  our  code  for 
extrapolation  to  the  limit  of  infinite  number  of  generations: 


*J<i  — o„+  | 


+  i  —  25„  +  S„ 


(10) 


where  S„s  are  partial  sums  of  the  image  series  and  is  an  improved  estimate  for  the  n-» oo 
limit. 

The  infinite  array  problem  is  equivalent  to  that  of  a  cylinder  placed  in  a  “frictionless  duct” 
with  walls  parallel  to  the  .x-axis,  of  width  equal  to  the  spacing  of  the  infinite  array.  In  this 
case,  an  initially  placed  doublet  at  the  center  of  the  cylinder  is  reflected  across  the  “duct 
walls”  and  the  image  generation  sequence  proceeds  as  in  the  case  of  the  finite  array.  We  had 
no  problem  in  converging  on  the  correct  solution,  that  contained  the  duct  walls  as 
streamlines,  for  S>3.  For  S<  3  however,  due  to  the  stronger  image  interactions  the 
streamlines  of  the  flow  field  generated  with  the  method  described,  did  not  exactly  coincide 
with  the  duct  walls  due  to  fluid  “leakage”  across  them.  To  correct  for  this  situation  we  found 
it  convenient  to  add  a  few  regularly  spaced  undetermined  singularities  on  the  y-axis.  The 
strength  of  these  singularities  was  determined  by  requiring  no  side  leakage  across  the  duct 
walls,  at  an  equal  number  of  equidistant  collocation  points,  placed  along  the  walls, 
extending  to  10R  from  the  origin. 

From  Levin’s  theorem  given  in  the  Appendix,  we  know  that  the  value  of  the  stagnation 
point  velocity  gradient  in  potential  flows  plays  a  crucial  role  in  determining  the  onset  of 
inertial  impaction  by  imposing  the  value  of  Stkc,|,,  below  which  no  capture  occurs.  It  is  then 
appropriate  at  this  point  to  examine  the  influence  of  Ncyl  and  S  on  the  dimensionless 
(normalized  with  Ux/R)  velocity  gradient  [du/d.x],p.  For  an  isolated  cylinder  the  result  is 
well-known:  [du/d.x]sp  s  =  3.  =  2.  ' 

Figure  2  depicts  the  influence  of  the  number  of  cylinders  in  the  array  on  [du/dx],p  of  the 
central  cylinder  for  a  fixed  spacing  5  =  4.  The  dashed  line  is  the  infinite  array  limit, 
[du/d.x]5p,  Nitl=r  =2.3534.  As  expected,  increasing  the  number  of  cylinders  creates  larger 
velocities  in  (he  vicinity  of  the  collector  due  to  the  smaller  area  available  to  the  flow  and 
?accordingly^ larger  gradients,  since  the  flow  must  decelerate  locally  from  a  larger  value  to 
zero  as  the  stagnation  point  is  reached.  Figure  3  shows  the  effect  of  S  on  [du/dx],p  for  the 
infinite  array.  This  figure  also  includes  experimental  results  (open  symbols)  extracted  from 
“blockage  factor”  measurements  oT  high-Re  momentum  and  heat  transfer  to  a  cylinder  in 
a  duct  given  in  Zukauskas  and  Ziugzda  (1985).  The  potential  flow  model  captures  correctly 
the  relative  trends  but  lies  consistently  below  the  measurements  as  if  the  experiments 
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Fig.  2.  Dimensionless  stagnation  point  velocity  gradient  of  central  cylinder  for  S  =  4  as  a  function  of 
the  number  of  cylinders  in  the  array.  The  dashed  line  is  the  infinite  array  limit. 


Fig.  3.  Dimensionless  stagnation  point  velocity  gradient  of  a  cylinder  in  an  infinite  array  as 
a  function  of  5.  Experimental  measurements  for  a  cylinder  in  a  duct,  extracted  from  Zukauskas  and 

Ziugzda  (1985). 


corresponded  to  a  higher  ‘’effective  blockage’’  (smaller  spacing)  than  actual.  A  possible 
cause  for  this  might  be  that  the  boundary  layers  (BLs)  on  the  cylinder  and  on  the  duct  wall 
increase  the  “effective  blockage”  by  displacing  the  adjacent  flow.fdri vin Athuj) the  stagnation 
point  gradients  higher.  Indeed^hen  we  correct  approximately  for  such  effects  by  subtract¬ 
ing  from  S  the  estimated  total  thickness  of  the  boundary  layers  on  the  cylinder  and  duct 
walls  (filled  symbols  in  Fig.  3)  the  agreement  between  measurements  and  theory  becomes 
much  better. 

Our  results  for  [du/dx]sp  for  the  central  (“lest")  cylinder  of  the  array  as  a  function  of 
spacing  are  shown  in  Fig.  4.  Interestingly  enough,  [du/dx],p  develops  a  maximum  for  S  near 
3.5  for  Ncy |  =  3  which  increases  in  magnitude  and  moves  toward  smaller  S  values  as 
Ncy |  grows.  The  non-monotonic  behavior  *een  in  Fig.  4  is  the  result  of  the  competition 
between  “flo\y  through”  (dominant  at  large  S)  and  “flow  around"  (dominant  at  small  S)  the 
arrays  containing  a  finite  number  of  cylinders.  This  effect  is;of  course^absent  for  the  infinite 
arrays,  where  “flow  around”  the  array  cannot  occur. 

The  non-monotonicity  in  [du/dx],p  for  finite  arrays,  implies  that  there  is  an  optimum 
spacing  for  the  capture  of  particles  gf  a  given  size  near  the  onset  of  inertial  impaction  on  the 
central  cylinder,  since  Stk  for  this  spacing  will  reach  its  critical  value  sooner  than  for  other 
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Fig.  4.  Stagnation  point  velocity  gradient  of  central  cylinder-iependence  on  S  for  various  arrays. 


spacings  (see;also;section  4).  This  fact  means  that  if  the  purpose  is  to  decrease  deposition  this 
optimal  spacing  must  be  avoided,  while  if  the  goal  is  to  maximize  collection  this  spacing 
should  be  selected.  The  graph  of  Fig.  4  also  suggests  that  is  possible  to  develop  a  continu¬ 
ously  adjustable  inertial  separation  scheme  using  a  cylinder  array,  the  spacing  of  which 
could  be  varied  by  a  mechanical  arrangement,  elTectively  changing  on-line  Stkciil  and 
consequently  the  collection  characteristics  of  the  array. 


3.  CAPTURE  EFFICIENCY 


We  employed  what  might  be  called  a  “finite-analytic”  technique  for  computationally 
efficient  integrations  of  the  particle  momentum  equation,  given  below,  because  many 
trajectory  integrations  had  to  be  performed: 


£ur  =  u-u„ 
Dt  ~  t  ' 


(il) 


The  essence  of  the  technique  is  to  use  a  locally-valid,  analytic  solution  for  particle  motion  in 
a  piece-wise  fashion,  so  as  to  be  able  to  take  large  time  steps  in  integrating  the  equation  of 
motion.  Our  code  employs  two  such  solutions,  appropriate  for: 

(1)  particle  motion  in  a  linear  flow  field  (see,a!so,  Fernandez  de  la  Mora,  1985);  and 

(2)  particle  motion  in  a  constant  flow  field  (see,e.g.,Fuchs,  1964). 

In  equation  (11)  we  use  Stokes'  drag  law  for  simplicity,  since  Israel  and  Rosner  (1983)  have 
already  shown  how  to  take  into  account  non-linear  drag  laws  in  Stkeff.  The  particles  are 
injected  in  equilibrium  with  the  flow  at  upstream  infinity: 

U.e,  (12) 

and  are  advanced  to  a  predetermined  distance  (x  =  -50R  is  an  acceptable  value  in  terms  of 
speed  of  execution  and  accuracy  in  the  calculation  of  capture  efficiency)  using  the  aforemen¬ 
tioned  analytic  solution  for  motion  in  an  (approximately)  constant  flow  field.  The  integra¬ 
tion  then  proceeds  using  the  solution  for  motion  in  a  piece-wise  linear  flow  field.  The 
capture  efficiency  is  evaluated  as  usual  using  a  limiting  trajectory  scheme  (see  e.g.  Fuchs, 
1964),  computed  iteratively  with  a  bisection  technique.  In  the  calculations  reported  here  we 
treat  the  particles  as  point  masses  and  neglect  the  phenomenon  of  interception,  an  assump¬ 
tion  certainly  valid  considering  the  relative  sizes  of  suspended  particles  and  tube  collectors 
occurring  in  the  heat  exchanger  "fouling  applications  of  interest  to  us  (dp/2R  <  10-4). 
Inclusion  of  the  interception  effect  is  straightforward  (e.g^Fuchs,  1964),  but  it  adds  an 
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Stk 


Fig.  5.  Comparison  of  present  capture  efficiency  calculation  to  the  results  of  Ingham  el  al.  (1989)  for 

S  =  5.35. 


Stk 

Fig.  6.  Capture  efficiency  of  central  cylinder  in  infinite  cylinder  arrays  as  a  function  of  Stk  at 

various  spacings. 


Fig.  7.  Capture  efficiency  of  central  cylinder  in  finite  and  infinite  cylinder  arrays  as  a  function  of  Stk 
at  various  spacings.  Starting  from  the  uppermost  solid  curve  S:  2.5.  3,  4,  5,  oo,  4,  3,  10,  2.5. 


additional  variable,  the  interception  parameter,  dp/£R)lo  the  parameter  space  of  the  present 
problem,  {Ncy S,  Stk}. 

In  Fig.  5,  we  compare  our  >;e,pfStkf  results  for  an  infinite  array  with  S  =  5.35  to  those  of 
Ingham  et  al.  (1989),  as  read  from  curve  B  in  Fig.  5  of  their  paper.  It  can  be  seen  that  the 
agreement  is  excellent,  and  justifies  the  asymptotic  approximations  that  Ingham  et  al.  (1989) 
employed  to  reduce  the  computational  burden  of  their  flow  field  calculation  [cf.  their 
equations  ($)-(!  2)]. 

We  summarize  our  results  for  infinite  arrays  at  various  values  of  S  from  2.5  to  oo  in  Fig.  6. 
As  expected  from  the  results  of  Fig.  4,  decreasing  S  makes  it  easier  for  particles  to  impact 
and  consequently  the  capture  efficiency  increases.  The  maximum  increases  occur  for 
Stk  ~(!?(1)  since  forStk-»oo  the  value  of  f/c,p  must  approach  1  irrespective  of  the  spacing  of 
the  array.  Figure  7  shows  the  capture  efficiency  of  both  finite  and  infinite  arrays.  Interesti- 
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ngly  enough,  the  finite  array  results  for  the  same  spacing  are  always  lower  than  the  infinite 
array  results.  Both  sets  of  results  follow  the  “ordering”  of  the  [dtt/d.x],p  curves  of  Fig.  4. 


4.  CORRELATIONS 

Motivated  by  Levin’s  theorem  (Appendix)  that  gives  Stkcril  in  terms  of  the  carrier  fluid 
stagnation  point  velocity  gradient,  we  define  an  effective  Stokes  number,  Stkcff  basing  the 
characteristic  flow  time  on  the  (dimensional)  stagnation  point  fluid  deceleration  rate, 
[du/d.v]-1: 


This  definition;in  combination  with  Levin’s  theorem,  maps  the  critical  Stokes  number 
(=xU *  /R)  for  each  array  configuration,  onto  the  isolated  collector  value  of  1/8.  Based  on 
these  considerations  we  conjecture  (see^also,  Israel  and  Rosner,  1 983)  that  this  definition  will 
also  make  possible  the  correlation  of  the  capture  characteristics  of  the  different  arrays  in 
terms  of  a  single  “universal”  curve. 

The  results  shown  in  Fig.  8,  provide  partial  support  for  this  conjecture.  The  open  data 
points  correspond  to  finite  regular  arrays  and  the  filled  points  to  infinite  arrays.  The  solid 
line  in  Fig.  8  corresponds  to  the  infinite  spacing  limit  (i.evthe  isolated  collector).  In  other 
words,  the  isolated  collector  »/clp  value  computed  at  the  prevailing  Stkefr  can  serve  as  an 
engineering  approximation  for  the  collection  efficiency  of  the  “test”  cylinder.  Comparing 
Fig.  8  to  Fig.  7  note  how  the  present  correlation  scheme  is  especially  successful  in 
compressing  the  data  in  the  vicinity  of  Stkcril  on  a  single  curve.  Although  the  spread 
increases  at  higher  Stkcff-values,  the  correlation  scheme  may  be  sufficient  for  engineering 
estimates  for  all  Stk-values  since  the  fractional  error  between  the  solid  curve  and  the  data  is 
decreasing  with  increasing  Stkcfr.  A  better  correlation  would  have  been  obtained  if  we 
eliminated  the  open  circles  which  correspond  to  a  closely-packed  (S  =  2.5)  three  cylinder 
array.  This  array  geometry  is  closer  to  that  of  a  ribbon  in  cross-flow  than  an  isolated 
cylinder,  and  the  flow  field  in  the  vicinity  of  the  test  collector  is  quite  different  from  the  other 
cases.  The  observed  spread  at  high  Stktff  is  due  to  the  fact  that  in  this  limit,  the  various  array 
flow  fields  experienced  by  the  particles,  are  more  dissimilar  than  the  “almost  universal” 
stagnation  flow  regions  experienced  by  slightly  supercritical  particles.  For  completeness,  we 
cite  the  Israel  and  Rosner  (1983)  correlation,  for  »;c,p-(Stkc(f}  (extended  here  to  the  case  of 
cylinder  arrays)  valid  for  Stkcff >0.I4t,  although;as  mentioned,  other  acceptable  alternat- 


Fig.  8.  Correlation  of  inertia!  impaction  capture  efficiency  of  a  central  cylinder  in  finite  (open 
symbols)  and  infinite  (filled  symbols)  arrays  in  terms  of  Stk,„.  Circles  correspond  to  S  =  2.5,  squares 
to  3,  triangles  to  4  and  stars  to  5.  The  solid  line  is  the  isolated  collector  curve. 


1  The  behavior  off)clp4Stk,„f  for  Stkcr(-*Stkc,h  is  non-analytic  and  is  discussed  in  a  set  of  unpublished  notes  by 
Prof.  J.  Fernandez  de  la  Mora  (1985). 
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Fig.  9.  Correlation  of  stagnation  point  impact  velocity,  [upJ].r  for  the  central  cylinder,  in  an  array 
of  ihrec  collectors,  at  5  =  2.5  (squares),  4  (circles),  oc-  (triangles),  in  terms  of  Stk<ff. 


Ftg.  10.  Correlation  of  stagnation  point  particle  concentration  enrichment  p-./p.,  for  the  central 

cylinder,  in  an  array  of  three  collectors,  at  S  =  2.5  (squares).  4  (circles),  oc  (triangles),  in  terms  of 

Slkcrr. 


ives  (Wang,  1986;  Wessel  and  Righi,  1988)  also  exist: 

,7c»p4Stkeff)-%f  I  +  1.25(Stkefr— 4)_  1  —  1.4  x  10 _ 2 (Stkerr  —  A)- 2 

+  0.508  xlO-4(Stktrf-i)-3]-’«.  (14) 

Correlation  of  the  data  in  terms  of  Stktfr  seems  to  be  possible  for  other  quantities  of 
interest  as  well  (at  least  for  an  array  of  three  collectors  tested  so  far), such  as  impact  velocities 
at  the  stagnation  point  (Fig.  9),  concentration  enrichment  factors  (Fig.  10),  etc.  which  will  be 
addressed  in  more  detail  in  the  future.  1 


5.  CONCLUDING  REMARKS 

Inertial  impaplion  calculations  for  linear,  regularly  spaced,  circular  cylinder  arrays  reveal 
that  the  principal  elfect  of  multiple  target  proximity  is  to  change  the  local  gas  deceleration 
rates  of  the  carrier  fluid,  thereby  shifting  Stkcril.  Thus,  our  numerical  results  for  particle 
capture  on  the  central  cylinder  in  arrays  containing  from  three  to  an  infinite  number  of 
flanking  cylinders,  are  fairly  well-correlated  using  the  relevant,  effective  Stokes  number, 
Stk,rr,  Le.  basing  the  “characteristic  flow  time”  on  the  actual  (  not  on  the  isolated-target) 
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gas  deceleration  rate  in  the  vicinity  of  the  stagnation  point.  Overall,  it  now  appears  possible 
to  make  engineering  estimates  for  the  capture  efficiency  of  fluid-dynamically  interacting 
cylinders  in  inviscid.  cross-flow  from  a  knowledge  of  the  stagnation  point  gradient,  given  in 
Fig.  4  as  a  function  of  number  of  cylinders  and  spacing,  and  the  well-known  correlations  for 
a  single  cylinder  (Israel  and  Rosner,  1983;  Wang,  1986;  Wessel  and  Righi,  1989).  Work  in 
progress  in  this  area  is  focusing  on  extending  the  correlation  approach  to  angular  distribu¬ 
tions  of  quantities  of  interest,  examining  other  cylinder  arrangements  (e.g.  staggered  arrays) 
as  well  as  studying  asymmetric  deposition  patterns  on  the  side  cylinders. 

It  would  also  be  of  interest  to  compute  (using  the  inviscid  velocity  distributions  on  the 
surface  of  the  collector)  the  fluid  mechanics  of  the  boundary  layer  (BL)  developing  along  the 
cylinder  surface,  to  estimate  collector  proximity  efTects  on  BL  mass,  momentum  and  heat 
transfer  rates  to  our  representative  cylinder  in  the-array.  We  note  however,  that  self-similar 
BL  solutions  based  on  the  inviscid  stagnation  point  velocity  gradient  (Fig.  4),  are  often 
acceptable  engineering  approximations  to  the  flow  up  to  ca  20°  from  the  stagnation  point  of 
the  bluff  collector  (Schlichting,  1955),  where  a  significant  part  of  deposition  also  occurs  in 
practice  (Rosner,  1986).  It  is  straightforward  then  to  estimate  deposition  rates  for  this 
geometry  under  the  simultaneous  action  of  inertial  (Stk <Stkcril)  and  thermophoretic/ 
diffusive  effects  applying  correlations  and  boundary  layer  deposition  theory  for  curved 
walls  (Konstandopoulos  and  Rosner,  1990;  Konstandopoulos,  1991). 

Furthermore,  while  we  have  deliberately  assumed  that  particles  arriving  at  the  surface  do 
not  rebound  upon  impact,  the  present  results  can  easily  be  corrected  for  effects  of  incom¬ 
plete  particle  capture  using  analytic  expressions  for  the  sticking  coefficient  of  particles 
impacting  clean  collectors  (Konstandopoulos,  1991;  Wang  and  John,  1988).  Along  these 
lines,  the  stage  is  now  also  set  for  extensions  of  the  calculations  performed  here  to  account 
for  hydrodynamic  interference  effects  due  to  thick  deposit  growth  on  the  collector 
(cf.  Rosner  and  Giines,  1993).  An  attractive  route  would  be  to  explore,  for  example,  conformal 
mapping  techniques  in  order  to  “deform”  the  non-circular  collector  cross-section,  resulting 
from  deposit  accumulation  on  the  cylinder,  back  to  a  circular  shape.  In  addition,  such 
calculations  can  now  incorporate  sticking  coefficients  of  particles  impacting  on  deposits, 
extracted  from  discrete-particle  simulations  of  deposit  growth  (Konstandopoulos,  1991; 
Konstandopoulos  and  Tassopoulos,  1990). 
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APPENDIX:  ON  THE  CRITICAL  STOKES  NUMBER 

It  has  been  know  n  for  quite  some  time  that,  under  certain  conditions,  there  exists  a  finite,  critical  Stokes  number, 
Stk„|,,  for  inertial  impaction  from  (lows  past  bluff  bodies.  For  Stk<Stkclll  no  particle  trajectory  which  is  in 
equilibrium  with  the  flow_at  infinity  reaches  the  surface  with  finite  velocity,  in  the  absence  of  other  transport 
mechanisms.  Taylor  [194QJ  obtained  the  value  of  Stk„i,  for  impaction  from  inviscid,  stagnation-point  flow,  making 
certain  assumptions  about  the  initial  conditions  of  the  particles,  while  general  conditions  for,  and  proof  of,  the 
existence  or  Stkcr{,  for  inviscid  flows  past  bluff  bodies,  were  summarized  in  an  important  yet  not  widely  known  (in 
•he  Western  literature)  monograph  by  Levin  (1961).  We  will  term  this  proof  Levin’s  theorem,  stated  below 
^paraphrasing  excerpts  from  Levin  (196l)jas: 

Levin  s  Theorem:  For  flows  past  axially  or  n-th  order  (n^  2)  rotationally  symmetrical  bodies  of  character¬ 
istic  dimension  L,  which  are  parallel  to  the  body's  axis  of  symmetry  at  infinity  (with  velocity  l/K),  and  the 
streamlines  of  which  terminate  at  the  upstream  stagnation  point  as  straight  lines  with  slope  a,  there  exists 
a  critical  Stokes  number,  Stkcth  =  rl/K/L=  l/(4a),  for  a  broad  range  of  initial  conditions,*  below  which  no 
particle  reaches  the  body. 

More  recently,  in  an  article  published  in  (his  Journal,  Ingham  ef  of.  (1990)  have  also  given  a  detailed  proof  or  the 
essence  of  Levin's  theorem  for  potential  flows  and  suggested  that  Slk„„  for  viscous  flows  is  zero.  They  reached  this 
conclusion  by  examining  particle  trajectories  in  the  quadratic  flow  u=bx 1  obtained  as  the  leading  term  of  the 
Taylor  expansion*  of  the  normal  component  of  the  fluid  velocity,  u  in  the  distance  x  from  a  boundary  surface. 
However,  Levin  ( 1 96 1 )  has  proven  the  existence  of  StkeI)1  >0  for  viscous  flows  as  well,  given  appropriate  initial 
conditions.  Levin  implicitly  recognized  that  in  this  case  the  onset  of  inertial  impaction  is  governed  by  phenomena 
occurring  in  a  region  where  the  quadratic  flow  is  a  poor  approximation  to  the  flow  field,  but  he  nevertheless 
suggested  a  scaling  relation:  Stkcrl,  =  const,  b  l,J,  and  evaluated  the  constant  by  comparison  to  numerical 
estimates  ofStkCI„  of  viscous  flows.  Although  Stkcr|,>0  has  been  reported  based  on  numerical  studies  for  some 
viscous  flow  cases  (Fuchs,  1964;  Levin,  1961),  it  is  also  possible  to  compute  Stkcf„  semi-analytically  following  the 
approach  of  Michlel  and  Norey  (1970).  These  authors  treated  carefolly  the  singularity  of  the  particle  momentum 
equation  at  the  stagnation  point  and  calculated  Stk„„  for  inertial  impaction  from  creeping  flow  past  a  sphere  to  be 
1-2 1 1 94.  Langmuir  [I948J  had  obtained  numerically  for  that  case  Stkclll%  1.214. 


*  Discussed  by  Levin  (1961);  see, also, Fernandez  de  la  Mora  (1980,  1985)  and  Ingham  et  al.  (1990). 

This  quadratic  dependence  is  a  consequence  of  the  no-slip  condition  of  the  steeMMaiM  component  of  the 
velocity  and  incompressibility. 


/  rr>  FORCE  OF  SCT'rrmc  P)-EF«r/.„ 

«;r  lini;:  ■;  nva  Jch 

h  "T  t 


(.AFSC) 


.......  f-u^lii^,  V <_  1 1; j 3 .9  p p D  inn  10 

/"M?  ^ tn  unlimited.  *  . 


and  is 

Deposition  of  particles  from  potential  (lows 


I  i l\;i  0  1 jJn  OQTaii]  jraritiaitQw  .  .... 

rrf  OTOrnpvfe  recall  that,  as  Fernandez  dc  la  Mote (19851  has  shown,  for  arbitrary  linear  and  locally  linear  flows 
characterized  by  a  deceleration  time  scale  to'  the  partieK  motion  becomes  critically  damped  by  the  host-flow 
when  the  product  tut  becomes  equal  to  l/[4niax(;.,)]  where  A,  stands  for  any  of  the  eigenvalues  or  the  velocity 
gradient  tensor  (rendered  dimensionless  with «.»).  Whether  the  Stk,  b  l/[4  max(/.,)]  thus  determined  coincides  with 
Stk,,i,  for  inertial  impaction,  depends  on  the  nature  of  the  flow  field.  The  two  numbers  coincide  for  potential  flows 
where  the  maximum  velocity  gradient  occurs  on  the  boundary  of  the  body. 

For  viscous  nows  however.  Stkt,„  can  be  significantly  higher  than  Stk.  due  to  the  fact  that  the  maximum 
velocity  gradient  occurs  away  from  the  collector  wall  and  the  particle  motion  has  the  chance  to  become 
overdamped  again,  inside  the  almost  stagnant  viscous  region  adjacent  to  the  wall.  For  example,  in  the  case  of 
creeping  flow  past  a  sphere,  the  particle  motion  becomes  first  underdamped  for  Stk,  >2/3  at  a  distance  (V2-  1)R 
from  the  sphere  surface  along  the  stagnation  line  where  the  maximum  velocity  gradient  occurs,  yet  Stkff„  for  this 
flow  as  we  saw*  above  is  1.21 194.  Similarly,  for  a  cylinder  in  Oseen  flow,  underdamping  at  first  occurs  when 
Stk,  >0.6495(2.002-  In  Re)'  at  a  distance  (N/5—  l)R  from  the  surface  along  the  stagnation  line.  Stk„„  for  this  case 
is  reported  to  be  4.3±0.l  for  Re  =  0.l  (Fuchs.  1964). 
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Use  of  a  Generalized  Stokes  Number  to  Determine 
the  Aerodynamic  Capture  Efficiency  of 
Non-Stokesian  Particles  from  a  Compressible 
Gas  Flow 

R.  Israel  and  D.  E.  Rosner 

High  Temperature  Chemical  Reaction  Engineering  Laboratory, -Chemical  Engineering  Department, 
Yale  University.  New  Haven.  CT  06520 


The  aerodynamic  capture  efficiency  of  small  but  non- 
diftuslug  particles  suspended  In  a  high-speed  stream 
filming  past  a  target  is  known  to  be  influenced  by 
parameters  gosernlng  (a)  small  particle  inertia,  (b)  de¬ 
partures  from  the  Stokes  drag  law  (associated  with  local 
particle  Reynolds  numbers  greater  than  unity),  and  (c) 
carrier  fluid  compressibility  (at  nonncgligible  free- 
stream  Mach  numbers),  fly  defining  an  effective  Stokes 
number  in  terms  of  the  actual  (prevailing)  particle  stop¬ 
ping  distance,  local  fluid  viscosity,  and  Invlscld  fluid 
velocity  gradient  at  the  target  nose,  we  show  that  these 


effects  are  well  correlated  in  terms  of  a  "standard” 
(cylindrical  cullcclor,  Slokes  drag,  incompressible  flow, 
Re,/2  :*•  1)  capture  efficiency  curve.  We  are  thus  led  to  a 
correlation  that  (a)  simplifies  aerosol  capture  calcula¬ 
tions  In  the  parameter  range  already  included  in  previous 
numerical  solutions,  (b)  allows  rational  engineering  pre¬ 
dictions  or  deposition  in  situations  not  previously  speci¬ 
fically  calculated,  (c)  should  facilitale  the  presentation  «f 
performance  data  for  gas  cleaning  equipment  and  aerosol 
instruments. 
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